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Abstract 

The study of the hydrodynamics of bubble growth in first-order phase transitions is 
very relevant for electroweak baryogenesis, as the baryon asymmetry depends sensi- 
tively on the bubble wall velocity, and also for predicting the size of the gravity wave 
signal resulting from bubble collisions, which depends on both the bubble wall velocity 
and the plasma fluid velocity. We perform such study in different bubble expansion 
regimes, namely deflagrations, detonations, hybrids (steady states) and runaway so- 
lutions (accelerating wall), without relying on a specific particle physics model. We 
compute the efficiency of the transfer of vacuum energy to the bubble wall and the 
plasma in all regimes. We clarify the condition determining the runaway regime and 
stress that in most models of strong first-order phase transitions this will modify ex- 
pectations for the gravity wave signal. Indeed, in this case, most of the kinetic energy 
is concentrated in the wall and almost no turbulent fluid motions are expected since 
the surrounding fluid is kept mostly at rest. 
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1 Introduction 

A cosmological first-order phase transition could have far-reaching consequences and could 
lead to many interesting phenomena, as for example electroweak baryogenesis [1J, primordial 
magnetic fields [2] or a stochastic background of gravitational waves [3H5]. All of them are 
based on the fact that a first-order phase transition proceeds by bubble nucleation and, in the 
process, a large portion of the available vacuum energy is stored close to the bubble walls. In 
all these phenomena, some essential characteristics of the first-order phase transition enter, 
namely the velocity of the expanding bubbles and the efficiency coefficients that keep track 
of the energy budget of the transition, that is, how the free energy of the Higgs field </> 
initially available is distributed among bulk fluid motion, thermal energy of the plasma and 
gradient/kinetic energy of the Higgs field. 

Quite generally, treatments of the bubble wall velocity assume that the expansion of the 
bubbles is hindered by some form of friction, such that the bubble wall reaches a constant 
speed after a rather short time of order ~ 1/m, where m represents the typical mass scale 
associated with the transition (i.e. electroweak mass scale in the case of the electroweak phase 
transition). Assuming that the free energy of the Higgs field is released into the plasma, a 
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hydro dynamic treatment of the plasma can be used to determine the fluid motion [HIE]- 
However, as long as the microscopic mechanism of friction is unknown, this approach leaves 
one free parameter, typically the wall velocity. In [7J it was argued that the velocity should 
be fixed by the Chapman- Jouguet condition, like in chemical combustion, and this condition 
was subsequently used in many studies of gravitational wave production j8]. The Chapman- 
Jouguet condition generally leads to supersonic wall velocities with a rarefaction wave behind 
the bubble wall. This scenario favors gravitational wave production, since fast moving walls 
are essential for the production of gravitational radiation in bubble collisions |i|[5|l9l4TT]. 
turbulence [121E13] or magnetic fields [T3] . 

Nevertheless, while being observed in chemical combustion, the Chapman- Jouguet condi- 
tion is unrealistic for cosmological phase transitions [H]. To replace the Chapman- Jouguet 
condition, friction has been studied through a phenomenological parametrization in the Higgs 
equation of motion [T5|[T6] . Extensive simulations of hydrodynamic equations in conjunction 
with an equation for the Higgs field have been used to identify the stable expansion modes 
of the bubbles [T ^fToTTIU] . It turns out that many solutions which are viable close to the 
phase transition front, where the free energy of the Higgs field is released into the plasma, 
are not stable globally. 

The friction coefficient appearing in the Higgs equation can be determined by solving 
Boltzmann-type equations close to the phase transition front [15j|20]. This approach leads 
to subsonic wall velocities in the Standard Model (SM) and its supersymmetric extension 
(MSSM) and is widely used in studies of MSSM electroweak baryogenesis [2TH25] . Subsonic 
wall velocities are crucial in electroweak baryogenesis, since this mechanism is based on 
diffusion of particle asymmetries into the plasma in front of the bubble wall and for a too 
fast wall, there is no time to build up a baryon asymmetry. 

In strong first-order phase transitions, the friction exerted by the plasma on the wall 
might not be sufficient to prevent the bubble wall from a runaway behavior in which the 
wall keeps accelerating, toward ultra-relativistic velocities, as pointed out recently in [27J. 
We also study this runaway regime and discuss the energy balance in that case. 

Our general goal in this paper is to present a unified picture of all the different regimes. 
We go beyond the (unjustified) Chapman- Jouguet assumption and provide general formula 
for the wall velocity, the fluid velocity and the efficiency factors accounting for the distribu- 
tion of energy among bulk fluid motion, thermal energy of the plasma and gradient /kinetic 
energy of the Higgs field. The paper is organized as follows. We review the hydrodynamic 
treatment of the plasma in sees. [2] and 13], and subsequently obtain the efficiency coefficients 
in sec. HI In sec. [5] we present a simplified treatment of the Higgs equation (a similar one 
has been carried out recently in We present detailed results both for the wall velocity 

and the efficiency coefficients that are required for the determination of gravitational wave 
spectra. We discuss the runaway regime in sec. |6] and summarize our results on the energy 
budget of first-order phase transitions in sec. [7J We conclude in sec. [8] and provide numerical 
fits to the efficiency coefficients in Appendix A. 
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2 Hydrodynamic relations 



In this section we introduce the basic concepts and set up the notation used for the hydro- 
dynamic analysis of the combined "wall-plasma" system [5J|7J[T5]. 

2.1 Basic concepts 

The energy-momentum tensor of the Higgs field <fi is given by 

"1 



2 d p <pdy-v (<p) 



where Vo(0) is the renormalized vacuum potential. The energy momentum-tensor of the 
plasma is given by 

^™ = E/(2^W*(M). (2) 

where the sum is carried out over the species in the plasma and fi(k,x) is the distribution 
function for each species. If the plasma is locally in equilibrium (perfect fluid) this can be 
parametrized as 

Tf— = wu,u v -g, vP) (3) 

where w and p are the plasma enthalpy and pressure, respectively. The quantity is the 
four-velocity field of the plasma, related to the three-velocity v by 

^ = -^= = (7,7v). (4) 

A constant background contributes to the total pressure [see eq. (1)] and from now on we 
will use p for this total pressure, including such contribution. 

The enthalpy w, the entropy density a and the energy density e are defined by 

where T is the plasma temperature. One then has 

w = e + p . (6) 

Conservation of energy-momentum is given by 

d^T^y = + ^T^ asma = . (7) 

We are primarily interested in a system with a constant wall velocity and, assuming there 
is no time- dependence, eq. (j7|) reads in the wall frame (with the wall and fluid velocities 
aligned in the z direction) 

d?T zz = d?T z0 = . (8) 
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Integrating these equations across the phase transition front and denoting the phases in front 
and behind the wall by subscripts + (symmetric phase) and — (broken phase) one gets the 
matching equations (in the wall frame): 

00 O O /\ 

w+v + 7 + + p + = W-V _7_ + p_ , w+v + 7 + = W-V _7_ . (9) 
From these equations we can obtain the relations 

P+-P- v + e_ + p + 

v+v- = , — = ■ . (10) 

e + — e_ V- e + + p_ 

To proceed further, one needs to assume a specific equation of state (EoS) for the plasma. 
2.2 Equation of state 

Usually the plasma is well described by a relativistic gas approximation. In the symmetric 
phase, 

P+=3 a + T +-e, e + = a + T*+e, (11) 

where e denotes the false-vacuum energy resulting from the Higgs potential (defined to be 
zero in the broken, T = 0, true-minimum phase). While in the broken phase 



P- = \a_Tt, e. = a.Tt, (12) 

with a different number of light degrees of freedom across the wall and therefore different 
values a + and a_ (with a + > a_) and different temperatures on both sides of the wall. These 
expressions correspond to the so-called bag equation of state. 

In the general case, the free energy (J 7 = —p, also called sometimes effective potential, 
including finite temperature corrections) of a plasma of particles with arbitrary masses mi(4>) 
is given, in the non-interacting-gas approximation, by 

d 3 p 



i 

where V is the T = effective potential, Ef = p 2 + mf , 

poo 

Y b/f (x)= dyy 2 \og[lTex P (-^x 2 + y 2 )], (14) 
Jo 

the — /+ signs in the last equation hold for bosons/fermions, the index % denotes the different 
species and Ni the internal degrees of freedom (that in this notation is negative for fermions). 
Species that are light compared to T behave as a relativistic gas while species much heavier 
than T are Boltzmann suppressed in the plasma and can be neglected. As is well known, 



4 



for small rrii/T the function Y b /j tends to a constant (— 7r 4 /45 for bosons and 77r 4 /360 for 
fermions) so that a± in the bag EoS are explicitly given by 

^ V .vs.- 7 

30 



a± 



light i 



(15) 



while for large rrii/T one has Yb/f(rrii/T) ~ ±exp (— rrii/T). Hence, only particles that have 
masses comparable to T can cause deviations from the bag EoS. In many cases this deviation 
is small because the free-energy is dominated by a large number of light degrees of freedom. 
For example in the Standard Model, the W ± and Z° bosons contribute to this deviation 
from the bag EoS (and to a lesser extent the tops) but altogether the modification in the 
energy density is at the percent level and hence small. 

If the deviations from the bag EoS are not so small one can still parametrize the plasma 
behavior in terms of quantities that mimic the bag EoS ones. In particular, from the free- 
energy (J 7 = —p) we can define 



a± = 



3 dp 



ATI dT 



3oj± 
47? 



e± = J(e±-3p±) 



(16) 



In the special case of bag EoS, these definitions reproduce the correct a± and give e + 
and e_ = 0. In terms of these quantities one can write 



P± 



-a±Ti - e± 



e± = a±Ti + e± 



(17) 



but now a± and e± are T-dependent quantities and should be interpreted with some care. 
By their definition (fTB"]) they satisfy: 



de. 



Tf da - 



&T± 3 dT ± 

There is a certain arbitrariness in the definition of a± but ours is especially convenient 
because some formulas we discuss in sect. [5] can be easily generalized even if the bag EoS is 
not valid. 

Using the bag equations of state (fTTj) and (fT2|) in eq. (fTD"]) we get 

1 - (1 -3a+)r 



where we defined 



v + v. 



3 
3 



3(1 + a H 
(1 — 3ct H 



l + 3(l + a+)r' 



a+Tf 



a_T 4 



(19) 



(20) 



The quantity a + is the ratio of the vacuum energy to the radiation energy density and 
typically characterizes the "strength" of the phase transition: the larger a + the stronger the 
phase transition. These two equations can be combined to give 



1 + a A 




6v- 



2 2 1 
+ a% + -a+ - - 



(21) 
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so that there are two branches of solutions, corresponding to the ± signs in eq. (j2ip . In the 
general case, with some deviation from the bag EoS, eqs. ffl9l) and fl2T]) still apply, with r 
defined as before and a + = (e + — e_)/(a + Tf). 

Figure [T] shows these solutions for several values of a + . The first branch (with positive 
sign) gives solutions of detonation type, v + > V-, corresponding to r < 1+ g ■ The function 
v + has a minimum at t>_ = c s . Notice that the point v + — t>_ = 1 is always a solution of 
this kind, with r = 1+ l a+ ; besides, in the limit r — > one obtains v + — 1 and t>_ = 1/3. 
The second branch of solutions exists only if a + < 1/3, and is of deflagration type, v+ < V- 
(it ranges over all values of vJ). Also for deflagrations, the function v + has an extremum at 
f_ = c s ; the endpoint, t>„ = 1, is given by the limit r — > oo and yields t> + = 3 ~^ + . The 
limiting case a + = corresponds to t> + t>_ = 1/3, connects both regions and will be relevant 
for the discussion of shock fronts in the next section. 

1 
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Figure 1: Contours of the fluid velocities v + and V- in the wall frame for fixed a+. In the shaded 
region in the top-left no consistent solutions to hydrodynamic equations exist. Flow profiles in the 
shaded region in the bottom-right decay into hybrid solutions, with V- = c~ (see sect. I3.3|> . 




To determine all relevant quantities in terms of the phase transition parameters, one still 
has to determine one of the two velocities, or equivalently, relate T + to T_; this will be the 
subject of sec. O Before doing so, we want to provide a more detailed description and a 
physical understanding of the above solutions as well as of the hybrid regime. Hence, for 
now, we treat the wall velocity as a free parameter, without any further assumptions about 
the microscopic properties of the plasma. In sections [3] and H] we present the analysis of 
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the different types of solutions for the motion of the plasma fluid. We will be particularly 
interested in the efficiency factor k (that measures how much of the vacuum energy goes into 
bulk kinetic energy and is a crucial quantity to determine gravitational wave production) as 
a function of the wall velocity. To define k, we first need the relativistic fluid equations for 
the plasma that we derive now. 

2.3 Relativistic fluid equations and bulk kinetic energy 

The starting point are the continuity equations (JTj) that, for a general plasma of form fl3]), 
read (see e.g. [6]) 

d^ v = u u d fl (u tl w)+u^wd ll u u -d u p. (22) 
Projecting along the flow then leads, using u^d^u^ = 0, to 

d^w) -u^p = 0. (23) 

Consider next the projection perpendicular to the flow with some space-like vector u = 
j(v,v/v) such that u^u^ = 0, u 2 = —1. This turns ([22]) into the relativistic Euler equation 

uV wd^u v - u u d u p = 0. (24) 

In the following we assume a spherically symmetric configuration. Additionally, because 
there is no characteristic distance scale in the problem, the solution should be a similarity 
solution that depends only on the combination £ = r/t where r is the distance from the 
center of the bubble and t is the time since nucleationQ £ is thus the velocity of a given 
point in the wave profile and the particles at the point described by £ in the wave profile 
move with velocity v (£), which is therefore the fluid velocity in the frame of the bubble 
center. This turns the different gradients into 

u li d» = -l(£-v)dt, « M 0" = !(l-£«)^. (25) 

Finally, this yields for eqs. (J2"3"]l and 

= 2^ + [l—fv{t-v)]dtv 



1 - vt — 



vO— = Ytt-v^v. (26) 



w 



The derivatives d^e and d(p can be related through the speed of sound in the plasma, 
c? s = (dp / dT) / (de / dT) , so as to get the central equation describing the velocity profile: 



2|=7 2 (l-<) 
with /i the Lorentz-transformed fluid velocity 



dp, (27) 



1 - £v 



1 The matching equations © can also be derived using energy-momentum conservation for similarity 
solutions in the rest frame. 
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In general, c 2 depends on the EoS for the plasma, being c 2 = 1/3 in the bag case. In the 
general case, c 2 will be ^-dependent, although in many cases of interest deviations from 1/3 
will be small. 

Eq. ( )27|) can then be solved (with the appropriate boundary conditions) to yield the 
velocity profile v{£) of the plasma. Subsequently, eqs. ( |26|) can be integrated to yield 



W {C) — w o ex P 



«(0 



1 + 



7 2 dv 



(29) 



In the calculation of the gravitational radiation produced in the phase transition one 
needs to compute the kinetic energy in the bulk motion of the plasma. We have now all 
ingredients necessary to perform such calculation. The ratio of that bulk kinetic energy over 
the vacuum energy gives the efficiency factor k as 



(30) 



where £ w is the velocity of the bubble wall. Notice that this definition coincides with the 
expression used in the gravitational wave literature, that is given by k = — ^3- J w v 2 ^ 2 R 2 dR, 
but differs from the definition used in ref. [5] by a factor 

We also numerically check energy conservation: Integration of Too over a region larger 
than the bubble (including the shock front) is constant in time, giving 



, 2 lx 3 
(7 - -)w - -w N 



3 

uj 1 



(31) 



where denotes the enthalpy at nucleation temperature far in front of the wall. This 
implies that the energy which is not transformed into kinetic bulk motion, but is used 
instead to increase the thermal energy, is 



1 - « 



3 



-(w - w N )£ 2 d£ 



3 



[e - e N )£ 2 d£. 



(32) 



3 Detonations, deflagrations and hybrids 

We can now use the previous fluid equations to describe the different kinds of solutions for 
the motion of the plasma disturbed by the moving phase transition wall. In the discussion 
below, the sound velocity in the plasma plays a very relevant role. This velocity will in general 
depend on £ and it is convenient to distinguish its asymptotic values in the symmetric and 
broken phases. We denote those two velocities by cf. In many cases, we expect the bag EoS 
to hold in the symmetric phase and therefore c+ = 1/ v^- 

Before embarking in the discussion of the different types of velocity profiles, it proves use- 
ful to study first in more detail the profile eq. f l2"T|) without worrying about physical boundary 
conditions. The different curves in Fig. [2] are obtained by solving for £ as a function of v 
[instead of the more physically meaningful v (£), the plasma velocity profile] using arbitrary 
boundary conditions and setting c s = l/y/3. This procedure has the advantage that £(i>) is 



a single-valued function. The meaning of the different regions will be explained later on. We 
see two different fixed points: one for £ = c s and v — and the other for £ = v — 1. This 
structure of fixed points can be understood analytically in a simple way by introducing an 
auxiliary quantity r to describe parametrically these curves in the plane (£, v). Using such 
parameter, eq. ( 1271) can be split in the two simpler equations 



^[(e-«) a -^(i-e«) 



(33) 



which clearly display the mentioned fixed points. The point (c s , 0) is in fact what is techni- 
cally called an improper node and all curves approach it tangentially to the v — line. One 
also sees that dv/dr > so that along each curve v grows monotonically with r. In the more 
general case with a ^-dependent c s deviating somewhat from l/\/3 one expects quantitative 
changes in these curves but the same qualitative behavior [with the fixed point at (cj",0)]. 




Figure 2: General profiles of the fluid velocity v(£) in the frame of the bubble center (with (? a = 
1/3). Detonation curves start below n(£,,v) = c s (dashed-doted curve) and end at (£, v) = (c s ,0). 
Deflagration curves start below v = £ and end at //(£, v)£ = c 2 s (dashed curve) corresponding to the 
shock front, as explained in the text. There are no consistent solutions in the shaded regions. 



Physical velocity profiles for expanding bubbles have to go to zero at some distance in 
front and behind the bubble wall and, in view of Fig. |2l this clearly requires some discon- 
tinuous jump, which usually can take place right at the phase transition front, where many 
other quantities also jump. The different types of possible velocity profiles that result are 
discussed next. 
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3.1 Detonations 



A pictorial representation of a typical detonation is depicted in Fig. [31 right plot. The 
corresponding velocity profile is as in Fig. HJ lower left plot. More precisely, in detonations 
the phase transition wall moves at supersonic speed £ w (£ w > c+) hitting fluid that is at rest 
in front of the wall. In the wall frame, the symmetric-phase fluid is moving into the wall at 
v + = £ w and entering the broken phase behind the wall where it slows down so that v _ < v + . 
In the rest frame of the bubble center, the fluid velocity right after the wall passes jumps to 
v (£,w) — ^{ v +i v -) (the Lorentz transformation (|28|) from the frame of the wall to the rest 
frame of the center of the bubble) and then slows down until it comes to a stop, at some 
£ < £ w , forming a rarefaction wave behind the wall. From the previous discussion we know 
that v will go to zero smoothly at £ = c~. 
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Figure 3: Pictorial representation of expanding bubbles of different types. The black circle is the 
phase interface (bubble wall). In green we show the region of non-zero fluid velocity. 



In order to obtain a consistent solution in the region c~ < £ < £ w , one needs < d^v < oo 
which, using eq. ff27j) . requires //(£) > fi(£ w ) > c~ behind the wall. Consequently, detonation 
solutions are confined to the lower right corner of fig. [2j as indicated. Boosting to the wall 
frame this implies v_ > cj, since f_ = v (£«,)). Therefore, detonations can be divided 
into Jouguet detonations (u_ = c~) and weak detonations (u_ > c^b strong detonations 
(v- < cj) are not consistent solutions of the fluid equations, see fig. do 

Fig. H] shows also the enthalpy profile (bottom right) for a detonation. Concerning this 
profile, remember that the matching conditions across the wall give 

where the subscript iV denotes the plasma at the temperature of nucleation far in front of 

2 As c~ can be different from l/v3 in the most general case, the forbidden region w_ < cj, shaded in 
Fig. []] will be shifted in those cases. 



10 



deflagration 



0.6 



0.4 



0.2 - 







"i 1 1 1 1 1 r 



cc + = 0.263 



r = 10.3 



0.45 0.5 0.55 0.6 0.65 



I 2 
q 

/l 



deflagration 
-i 1 1 1 1 1 1 1 r 



w w 






0.45 0.5 0.55 0.6 0.65 



0.6 



0.4 



0.2 



0.6 



hybrid 



1 


1 


oc + = 0.052 




r =2.06 













0.65 



0.7 



0.75 



D 
-a 

d 2 



~% 1 



0.6 



hybrid 



w w 



\ sh N 



0.65 



0.7 



0.75 



detonation 



detonation 



> 



0.5 
0.4 
0.3 
0.2 - 
0.1 - 




0.5 



a + = 0.091 
r = 0.458 




0.8 




Figure 4: Examples of the fluid velocity (in the plasma rest frame), enthalpy and entropy profiles 
for a subsonic deflagration, a deflagration with rarefaction wave (hybrid) and a detonation, for 
a_/a + = 0.85. The bubble of broken phase is in gray. For detonations, the fluid kinetic energy 
and thermal energy are concentrated near the wall but behind it i.e. inside the bubble, while they 
are located outside (mostly outside) of the bubble for deflagrations (hybrids). 



the wall. Then, eq. ( )29l) transforms into 



w{i) = w N 



1 

Si, 



exp 
11 



7 2 n dv 



(35) 



Similar formulas can be derived for other quantities like the entropy (also shown in fig. HJ), 
the temperature, etc. It is straightforward to show that 

T_>T + =T N . (36) 

by using the detonation condition r < 1/(1 + 3aq_). 

3.2 Deflagrations 

A pictorial representation of a typical deflagration is depicted in Fig. [3J left plot. The 
corresponding velocity profile is as in Fig. HI upper left plot. In contrast with detonations, 
in deflagrations the plasma is at rest right behind the wall, so that the wall velocity is now 
£ w = u_. These solutions correspond to the lower branches in Fig. [TJ with a fluid velocity 
that is larger behind the wall than in front, V- > v + . From Fig. [I] we also see that in this 
case the hydrodynamic relations across the wall imply v + < c s . The fluid velocity just in 
front of the wall jumps to v(£ w ) = fi(v-,v + ). As t>_ > v + , we get v(£ w ) < £ w , so that the 
profile of deflagration solutions start below the line v = £, as indicated in fig. [2J 

As one moves out in £, decreases and eventually would become double- valued before 
reaching zerclfl. Now we cannot accommodate the required jump to zero velocity using the 
phase transition front, which is already fixed at the beginning of the velocity profile. The 
way out of this dilemma is that the flow can drop to zero in a shock- front: we can use eq. (1191) 
where now quantities with subscript ± denote the two sides of the shock-front instead of the 
phase transition wall. (Notice also that the shock occurs in the symmetric phase where in 
general we can assume the bag EoS holds.) At the shock-front there is no discontinuity in 
e, and therefore a + = 0, yet a discontinuity in velocities is possible. As given by eq. ( IT9|) 
and shown in Fig. [TJ this requires [17] that the inward and outward velocities of the fluid 
in the shock-front rest-frame fulfill the relation t> + t>_ = 1/3, which in the plasma rest frame 
translates to fi(£, s h,Vsh) £,sh = c 2 s . This condition determines the position of the shock front, 
which can be seen to occur before the singular point = c s (where dv/d^ — > — oo) 

is reached, see fig. [2j 

Fig. HJalso shows the enthalpy and entropy profiles (top right) for a deflagration. Notice 
that the temperature T + in this case does not coincide with the temperature of the plasma 
outside the shock T/v, since thermodynamic quantities are changing in front of the wall and 
also at the shock front discontinuity. [For the same reason, explicit formulas like (j3"5|) would 
be now more involved.] As illustrated by the plots in Fig. H] one always has the inequality 
u + > un, i-e. a + T^ > a^Tfj. In fact, for deflagrations one always has the inequalities 

T + > T sh > T N > T_ (37) 

although the last one depends on the details of the Higgs equation of motion, discussed in 
sect. [5j Hence the limit a + < 1/3, that one found for a + = e/(a + T+) in deflagrations, 
translates into a weaker bound on the ratio = e/(a N T^) as > a + . 

An upper limit on will be important later on to set an upper limit on 77. How large 
can «tv be for a fixed wall velocity? From a + < 1/3 we get the upper bound < w + /(3a;jv) 

3 In spite of appearances, this happens even for very low values of the starting v(£) due to the improper- 
node nature of the point (c s , 0) (see discussion above). 
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and the ratio cannot be made arbitrarily large. Inspection of the deflagration profiles 

in Fig. [2] shows that will be maximized by the case with strongest shock-front, which 

corresponds to the highest possible v (£«,), that is, v(£ w ) — Cw 111 the wall frame this gives 
v+ — and, using the matching conditions (El), this case also has U- — > 0, i.e. aJFt — > 0. 
Physically this would represent a limiting case for which the transition is such that the broken 
phase inside the bubble is empty: the plasma is swept away by the wall (thus leading to the 
strongest possible shock-front) and larger values of rj cannot be realized microscopically. The 
same matching conditions also give us u + — > 4e, which represents a fixed upper bound for uo + . 
Going now to the shock rest frame and matching there one gets un — cu^/jC^l — £,sh)/(^sh — ^s)~ 
The quantities u s h and ^ can be obtained once the boundary condition v (£„,) = £ w is fixed 
and lead to the minimal value of ujn- The resulting upper bound on can be fitted 
numerically as a function of the wall velocity, and one gets 



max I ^ 
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as derived in Appendix |A] 
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Figure 5: Fluid velocity profiles (in the plasma rest frame) for deflagrations, hybrids and detona- 
tions, for different wall velocities and Ojv = 0.3. 



3.3 Hybrids 

From the previous discussion of velocity jumps it is clear that it should be possible to com- 
bine detonation and deflagration solutions into a new velocity profile that is a superposition 
of both types (thus the name hybrid), provided the wall is supersonic. In fact, it is known 
from hydrodynamics simulations [19] that supersonic deflagrations are not stable but develop 
a rarefaction wave, identical to the detonation profile discussed earlier. A pictorial repre- 
sentation of a typical hybrid bubble is depicted in Fig. El central plot. The corresponding 
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velocity profile is shown in Fig. HJ middle left plot. Applying the hydrodynamic constrains 
for deflagration and detonation solutions we have (see fig. [I]) fi(£ w , v(£~)) = v _ > c~ and 
l^{£,w, v (£,w)) = v + — c t ■ On the other hand, entropy considerations enforce t>_ < c~, so 
that the rarefaction wave has to be of Jouguet type, with t>_ = c~ . 

Therefore in these hybrid solutions the wall velocity £ w is not identified with either v + 
or and the phase transition front is followed by a rarefaction wave of Jouguet type 
(6uj > v - — C <T > v +) an d is preceded by a shock front. Typical enthalpy and entropy 
profiles are shown in the middle right plot of Fig. HI As the wall velocity increases, the 
deflagration part of the solution becomes thinner, and eventually disappears, as shown in 
Fig. El 




Figure 6: The flow velocities in the wall frame, v+ and V-, the velocity of the shock front £ s /j and 
the ratio o + /ojv as a function of the wall velocity (for aw = 0.1). 



In summary, one obtains the following picture. For subsonic wall velocities, the bubble 
expansion proceeds in general by deflagrations. If for constant oln the wall velocity is in- 
creased beyond the velocity of sound, the flow profile develops a rarefaction wave. If the wall 
velocity is further increased, the shock becomes thinner until it completely vanishes and the 
bubble expansion proceeds by a Jouguet detonation. When the shock front vanishes, some 
of the quantities right in front of the wall experience a jump, namely the temperature T + , 
the flow velocity v+ and also a + as illustrated in Fig. [6] Beyond this point the solutions are 
given by weak detonations (v- > c~). The hybrid solutions then fill the gap of wall velocities 
between Jouguet deflagrations and Jouguet detonations [19], as seen in Fig. [7] which gives 
the maximal flow velocity in the plasma rest frame as a function of the wall velocity £ w for 
fixed values of a^. For deflagrations and hybrids the maximal flow is in front of the wall 
while for detonations it is behind it. Notice that, in the transition from hybrids to detona- 
tions, the maximal flow velocity jumps without a gap in £ w . The whole flow profile can be 
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reconstructed using Fig. [2j 



4 Efficiency coefficients 

In the following we present the analysis of the efficiency coefficient k (the ratio of bulk kinetic 
energy to vacuum energy) following ref. [6l[lTJ,[5] . Steinhardt argued in his seminal work [7] 
that the relative plasma velocity behind the wall should just be given by the speed of sound, 
V- = c~ . Using this piece of information in eq. (1191) gives immediately 



\/a + (2 + 3a + ) + 1 
£,w = O = — 7= (Jouguet detonations), (39) 



\/3(l + n 

and 



k ~ —. = (Jouguet detonations). (40) 

0.135 + y/0.98 + a N v 1 y ' 

Notice that this differs from the result in ref. [5] due to our different definition of the efficiency 
factor k [see discussion after eq. f l30|) ]. 

As already mentioned, the Jouguet condition turns out to be unrealistic in cosmological 
phase transitions [19], and therefore one needs to rederive a formula for k that supersedes 
eq. ( HU1) . In the next section we will discuss how to relate the two velocities in the plasma 
v + and f_ (which is equivalent to the determination of the bubble wall). Knowing the wall 
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Figure 8: The efficiency coefficient k as a function of the wall velocity £ w for fixed a at. The 
dashed and dashed-dotted lines mark the transitions from deflagrations to hybrids and further to 
detonations. The dashed-dotted line corresponds to Jouguet detonations (the only case used in 
the literature, although with a missing factor). Analytical fits for K(aN,£,w) are provided in 
Appendix [S] 



velocity £ w and the parameter a + (or a^) the velocity profile is determined and k can be 
calculated using eq. (1501 independently from further assumptions on friction and microscopic 
physics in the plasma close to the wall (which are relevant to fix £ w ). 

The results are shown in Fig. [8] which gives k as a function of the wall velocity for 
several values of the vacuum energy a at. Note also that for large values of a^, small wall 
velocities are impossible, see Fig. [7] and the discussion about deflagrations in Sec. El The 
efficiency increases with ct^ and is maximal for the hybrid solutions. Nevertheless, according 
to numerical simulations, the detonation solutions are the only supersonic modes that are 
globally stable for small values of and realistically the maximal efficiency corresponds 
to the Jouguet case in this regime. The gravity wave literature focused on the Jouguet 
detonations (dashed-dotted line) and hence overestimated the efficiency k. However, we 
stress that this effect is mostly compensated by the missing factor £^ in the formula of n(a) 
we mentioned before. In appendix [Aj we give fits to the efficiency k shown in Fig. [8] as a 
function of the parameters and £ w . 

Finally, it is interesting to estimate the thickness of the plasma shell near the bubble 
wall where the kinetic energy in the plasma is concentrated (as this is relevant for GW 
production). In Fig. |H]we give the thickness A£ of a shell around the transition wall such 
that it contains a given fraction of the kinetic energy, Ak/k, as indicated. For each type 
of bubble we choose the shell so as to maximize that fraction. That is, for detonations the 
shell is (£ w ,£w — A£ w ); for deflagrations (£ w — A£ w ,£ w ) and in-between for hybrids. The 
dots in each line mark the boundaries between different regimes (deflagrations, hybrids and 
detonations, in order of increasing wall velocity). We see that, especially for hybrid solutions, 
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Figure 9: Curves with a fixed fraction of kinetic energy Ak/k in a shell of thickness A£ for 
atN = 0.1 and a at = 1 as a function of the wall velocity £ w . The lower black (upper brown) lines 
denote Ak/k = 0.5 {Ak/k = 0.9). The dots in each line mark the boundaries between different 
regimes (deflagrations, hybrids and detonations, in order of increasing wall velocity). 



the kinetic energy is localized in the shock, which can be rather thin in this case. For weaker 
deflagrations, the shock is thicker. For detonations, part of the kinetic energy is always in 
the tail of the rarefaction wave. 

In the limit of ultra-relativistic wall velocities (£ w — > 1) the velocities v± approach 1 with 
a relative slope v + /v_ that depends on a + (see Fig. [T]). Generically, the fluid velocity will 
not be ultra-relativistic in the plasma frame and the bubbles will expand as detonations (see 
Fig. [2]). Expanding f l2T|) in powers of l/y± one gets 

-t± = l + 3a N , v(£ w ) = fi(v+,v-) = 3aiV , (41) 
71 2 + 3a N 

where we have used a + = as we are in the regime of detonations. This relativistic case 
requires r = 1/(1 + 3a + ), which translates into p + — p_ = 2e. Numerically, the resulting 
efficiency in this limit is given by 

k ~ — ■ (relativistic £ w ). (42) 

0.73 + 0.083^0^ + a N v s 7 v ' 

This result is relevant for the analysis of the case of runaway walls in section [6j 



5 Bubble wall velocity for steady state walls 

When a bubble of the broken phase is nucleated, and is large enough to start growing, it 
will expand in an accelerated way, with the difference in free-energy across its wall acting as 
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driving force. There is however a resistance to this expansion from the surrounding plasma, 
which exerts a friction force that grows with the velocity of the moving wall. Eventually, an 
equilibrium between these two forces is reached after a short time of expansion and, since 
then on, the bubble wall keeps expanding in a steady state at a constant terminal velocity. 
As explained in the last sections, hydrodynamics alone cannot be used to determine this 
terminal wall velocity and one has to analyze the mechanism of entropy production and 
friction in the wall. 



5.1 EoM for the Higgs field and the friction parameter 77 

We take into account entropy production and friction through the equation of motion of the 
Higgs field 

^S + ??/pfe*»-°' (43) 

By decomposing 

ft(p) = f?(p) + 8fi(p), (44) 

where f^ q = l/[exp (Ei/T) =F 1] is the equilibrium distribution function of particle species 
i with Ef = p 2 + m 2 , eq. (j4"3l) takes the simple form (see also ref. [16J and more recently 
ref. [26]) 

{ty + — - JC{4>) = , (45) 
where the second term gives the force driving the wall and /C(0) stands for the friction term 

Friction is therefore due to deviations of particle distributions from equilibrium. In prin- 
ciple, calculation of JC((f>) requires solving a coupled system involving Boltzmann equations 
for particle species with a large coupling to the Higgs field. This intricate calculation has 
been performed in the Standard Model [TB] and in the MSSM [2U] and under the assumption 
that the deviation from thermal equilibrium is small, i.e. Sfi(p) <C fi(p), which is only true 
for weakly first-order phase transitions. 

In this paper, we want to follow a more phenomenological and model-independent ap- 
proach. In refs. [HH^B] a particularly simple choice for /C(</>) was used: 

K(<p) = T N fju' l d lt (f>, (47) 

(where Tjv is inserted just to make fj dimensionless) . This Lorentz invariant choice is mo- 
tivated by similar approaches in the inflationary context but, as we will see in the next 
section, it does not lead to the correct behavior for highly relativistic bubble wall velocities: 
this friction force could increase without bounds, due to the 7 factor appearing through 
u^d^cj), but we know from ref. [27] that at large wall velocities the friction term approaches 
a constant (see next section). 

Friction comes from out-of-equilibrium effects and the assumption that it depends lo- 
cally only on the plasma four-vector and a Lorentz scalar rj is too simplistic. In our 
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phenomenological approach we ensure that the friction force grows with v and not ^v. Such 
behavior could arise from a friction term in the Higgs equation of motion of the form 

(48) 

where the Higgs background is parametrized by a four-vector A M [such that (f){\^x^) and 
is (0, 0, 0, 1) in the wall frame]. One can show that the entropy production from such a term 
is always positive, as it should be. 

Assuming then that in the steady state the bubble is large enough so that we can use 
the planar limit, using (l4~8j) in eq. (I45p we get, in the wall frame, 

dT 

d 2 z <p- — = -T N fjvd z <P, (49) 

where z is the direction of the wall velocity. Note that the right-hand side would be multiplied 
by 7 if we use f )4T|) instead of f l4"g|) . If we multiply this differential equation by d z 4> on both 
sides and integrate across the wall, we get 

F)T r 

dzd z <p — = T N fj / dzv {d z <pf . (50) 



The integration of the force term could be simply performed if the free energy J 7 did not have 
an implicit dependence on z via the change in the temperature, T(z), with T(±oo) = T±. 
Using dJ-'/dz = (dJ 7 / d(p)d z (p + (9J 7 '/dT)d z T, one can rewrite the driving force of the bubble 
expansion as: 

F dr = J dzd z <j> ^ = Jt - J dzd z T^, (51) 

and, using e± and a(z) as defined in eqs. (|T6l) and ( TPT|) . one gets, without making assumptions 
on the plasma equation of state: 



7 dr = e + - e_ - - 



daT 4 . (52) 



By making further use of the definition of a(z) and assuming that the distribution functions 
for particle species are the equilibrium ones one can rewrite eq. fl52|) as 



where A Vo is the T = part of e + — e_ , that is, the difference in (T = 0) potential energy 
between the symmetric and broken minima (e, for the bag equation of state). This expression 
for the driving force will be useful in sect. EJ 

Notice that this force does not coincide with the latent heat A, given by 

A = e+ - e _ = ( e + a T 4 )|+ . (54) 
nor with the free energy (pressure) difference 



AJ 7 = v- ~ P+ = ( e - g T 



a „ 4 



(55) 
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Using the fact that the number of effective degrees of freedom decreases continuously in the 
wall, the last integral in ( 152]) is bounded as 



3 



If 1 

[a + -a_)T^ in < - J daT 4 < -(a+ - aJ)T^ ax , (56) 



where T min = Min{T_,T + } and T max = Max{T_,T + }. For weak phase transitions one 
typically has T + ps T_ and the concrete choice of the profiles in the wall are not important. 
For very strong phase transitions, the free energy is dominated by the vacuum contribution 
and the plasma contribution is rather small altogether: 

rp4. 

F dr ps AT ps e + — e_ — (a + — a_) — (for weak phase transitions) (57) 
F dr ~ AJ 7 PS A ps e + — e_ (for strong phase transitions) (58) 

For simplicity, let us approximate the integral in eq. (|56|) by the expression involving T + . 
For other approximations the results would not change qualitatively. In this case, the Higgs 
equation of motion gives 

^-k{ i -^hW' bv [d -" ? ■ (59) 

In the SM, the second term in (159]) is roughly (1 — a_/a + ) /3 = 0.05. For small values of a + 
the left hand side of ([59"]) is negative and no bubbles can nucleate (for nucleation one can 
consider T + = T_, in which case Fd r = J~+ — J 7 -, so that the wrong sign of the free-energy 
difference is the reason that prevents the phase transition). For larger values of a+, the 
left-hand side of the equation is positive and has to be balanced by the friction force so as 
to obtain a constant wall velocity, as already explained. We will rewrite the right hand side 
of ([59D as 

! dzvidAf^V—iv) , (60) 

which serves as the definition of 77. Here (v) denotes the fluid velocity average across the 
wall (in the wall frame), that we approximate as 

f dz v (d z d>) 2 1 . . . 

In this way, (!59]) is simply written a^ 

«+ - Ul ~ — ) = V— (v) ■ (62) 



3 \ a + J OtN 

Eq. (162]) reproduces the correct behavior for small and very large wall velocities, even 
though the parameter 77 will not be the same in these two limits. This simple phenomenolog- 
ical approach already reproduces almost all qualitative features found in the hydrodynamic 



4 Compared to ref. |26j . our velocity average does not include a 7 factor and we use a different expression 
for the driving force. 
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simulations performed in ref. [TO]. Notice that the right hand side of (I59p is an increasing 
function of the wall velocity (since v + and t>_ are, see also Fig. [6]). On the other hand, the 
left hand side of the equation is decreasing for small velocities and jumps when the flow 
solutions change from deflagrations to detonations (also the velocity v+ jumps but this is 
less essential). Hence, for a certain range of values of the friction coefficient 77, eq. (}59|) has 
two solutions while in ref. [16J three solutions were found. The discrepancy is due to the 
fact that we neglect the thickness of the wall compared to the shock and this would smooth 
out the jumps in v + and a + /a^, thus producing a third solution. The three solutions would 
in the present model then be a deflagration (or a hybrid solution), a weak detonation and a 
configuration close to a Jouguet detonation, in agreement with the results found in ref. [16J. 




Finally, at the critical temperature at which the two phases have degenerate free energy, 
one expects that the gravitational wave signal vanishes. However, an analysis using the 
Jouguet condition still implies supersonic bubble expansion. This led to the fact that the 
difference in free energy is often confused with the vacuum energy e in the literature on 
gravitational wave production that is based on the Jouguet condition. Replacing the Jouguet 
condition by the equation of motion (EoM) of the Higgs field solves this problem consistently, 
since at the critical temperature one obtains = (1 — a_/a + )/3 and, according to equation 
f l62]) . the bubble cannot expand. 

5.2 Wall velocity in the (77, a/v) plane 

A fully numerical calculation of (v) in a given model would require the following proce- 
dure. The starting point is the computation of the free-energy (in some approximation) as 
a function of the Higgs field and the temperature, J-"(0, T) = —p. First, the nucleation tem- 
perature T/v of the phase transition should be determined (this is standard). To obtain the 
steady-state profiles across the phase-transition wall in the planar limit of quantities like the 
velocity, temperature and Higgs field, one should integrate the following system of coupled 
differential equations: 



In addition to the first Higgs equation which we already discussed, the two extra equations 
correspond to the differential (and static) form of energy-momentum conservation. Their 
integration across the wall gives immediately the Steinhardt's matching conditions (j9]). 

One wants to solve the system (1631) for the primary quantities (j>(z), T(z) and v(z). The 
boundary conditions for <p(z) are 



where </>o(T!_) is the Higgs vacuum expectation value in the broken phase at some temperature 
T_ to be determined. The z = ±00 boundaries correspond to the assumption that the width 
of the phase transition front is much smaller than the width of the shell with non-zero 




%<f> ~-qT + T N f]vd z ^ = , 
d z [u>y 2 v] = , 

d z ^(^f + W7 V+p =0. 



(63) 



<p(— 00) = O (T_) and 4>(oo) = 



(64) 
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Figure 10: Contour plots of n and £ w as functions of r\ and (for a_/a+ = 0.85). The blue lines 
mark the transition to regions without solutions. The green lines mark the boundaries between 
stationary and runaway solutions. The red lines mark the transition from subsonic to supersonic 
deflagrations (hybrids). We superimposed the detonation region in the lower plots as a gray band. 



plasma velocity, which in general is a very good approximation. For fj fixed, the boundary 
conditions (say at z — — oo) for T(z) and v(z) cannot be chosen freely: e.g. if one fixes 
T(+oo) = T + (in general different from T N ) only one particular t>(+oo) = v + is selected 
and then all profiles <f)(z), T(z), v(z) can be determined. Detonation solutions will have 
v(+oc) — v+ — £ w > v(— oo) = t>_ and one should choose T(+oo) = T N . Deflagrations 
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will have instead v + < v _ = £ w and T(+oo) = T + has to be set consistently with such wall 
velocity (and corresponding deflagration hydrodynamic profile). Finally, hybrid solutions 
have v + < v _ = c~. Their wall velocity, £ w > c~ is not fixed unequivocally by v± but 
should match the choice of T+. Our analytical formula ( |62l) is expected to give a reasonable 
approximation to the more complicated numerical approach just described. Such simple 
formula is very useful to investigate in a model-independent way the parametric dependences 
of the wall velocity (and derived quantities) in different regimes of bubble expansion, which 
we do next. 

Armed with eq. ( 1621) we can calculate the wall velocity for given values of a a? (that 
measures the strength of the transition) and 77 (that measures the friction of the plasma). 
The results are shown in Fig. UHl which plots contour lines of the wall velocity £ w and the 
efficiency coefficient k in the plane (77, ajv), for the particular case of a_/a + = 0.85 (as in 
the SM). Let us first explain the different boundaries of the parameter space available in 
the plane (77, ajy). For the bubbles to be able to grow, the left hand side of fl62|) should 
be positive. This requires a + > (1 — a_/a + )/3 which, in this particular example, gives 
&n > a + > 0.05. Both the low horizontal boundary (a^r > 0.05) and the boundary at larger 
«jv and small 77 come from this requirement. That bubbles might not be able to grow even 
if the transition is in principle quite strong (large a^) is due to the fact that one can have 
a + C ajv for hybrid solutions (see figs. 161 and [T4|) . The boundary at > 1/3 and large 
77 corresponds to the extreme case with v + — > and T_ — > and determines the maximal 
possible value of the friction coefficient for a given ajy. We have plotted separately the cases 
of detonations and deflagrations. Runaway solutions are realized for strong transitions and 
small friction as one would expect. We see that for weak phase transitions and small friction, 
deflagrations compete with detonations (superimposed in the lower plots as a gray band). 
In the overlapping region, the wall velocity for the deflagration solution is always smaller 
than for the detonation solution, which could have indicated that the deflagration solution 
is reached first and is therefore realized. However, hydrodynamic simulations indicate that 
these solutions turn out to be unstable globally and in this regime the phase transition 
proceeds by detonation bubbles [19] (or eventually by runaway behavior if the upper bound 
on friction would be incorporated in the model). Henceforth, the only way of realizing 
supersonic deflagrations is in a regime where hybrids but no detonations are possible [19] . 



5.3 Microscopic determination of 7/ 

In the following we provide the connection of our results to the work [15j[20] in which the 
friction in the wall was explicitly calculated in the SM and MSSM@. The friction force in 
ref. [20J is of the form 

F fr = fj(v)T N j dz<f) 2 (<9 Z 0) 2 , (65) 

where 77 is a constant that depends mostly on the particle content of the model and its 
couplings to the Higgs. <f)(z) = \4>n\\ + tanh.(z/l w )] where l w is the wall thickness and (p^ 
is the Higgs vev at nucleation temperature. Comparison with our friction force shows that 

5 We consider those results more reliable than the ones of refs. [28ll26j which deduce particle distribution 
functions Sf without taking interactions into account. 
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our parameter r\ can be written as 

10 a + T N l w \T N J 

The coefficient fj was determined in the SM [15] (f) ~ 3) and in the MSSM [20] (fj < 100 with 
a sizable dependence on tan/3). A particularly interesting case is given by the parameter 
region of the MSSM that allows for viable electroweak baryogenesis. The bound on sphaleron 
wash-out implies 4>n/Tn > 1 and using TnI w 10, fj 100 one finds rj ~ 1/30. Due to a 
small difference in free energies, this leads to subsonic wall velocities (v ) = 0.05 ^0.1 [20] 
as required for the diffusion of CP-violating particle densities into the symmetric phase in 
front of the wall. This corresponds to a very weak phase transition with a value of just 
slightly above its lower bound (that depends on a_/a + ). Note that for models with a similar 
particle content the friction rj is not expected to change much, while the strength of the 
phase transition can increase significantly. This is for example the case in singlet extensions 
of the SM and MSSM which can easily lead to detonations or runaway solutions. 

In this section we have assumed that the bubble wall reaches at some (not too late) stage 
of the expansion a constant velocity. In this case the fraction of energy transformed into 
kinetic energy of the Higgs field becomes negligible, since it only scales with the surface of 
the bubble, while the similarity solutions of bulk motion scale with the volume. This can 
change in cases in which the wall keeps accelerating without reaching a terminal velocity, as 
discussed in the next section. 



6 Runaway walls 

It was recently argued [27J that the friction exerted on the Higgs wall by the plasma might 
be too small and the wall might continuously accelerate. In this case a constant fraction of 
the free vacuum energy is transformed into kinetic and gradient energies of the wall. In this 
section we analyze the energy balance and the efficiency coefficient in this situation. 

Let us first quickly present the main result of j2Z] that is based on the analysis of refs. [29- 
[31]. The passing phase-transition wall disturbs the distribution functions of particles in the 
plasma. As discussed in the previous section, if we knew such non-equilibrium distributions, 
fi(p,z), for each particle species, we could write, for the total force acting on the wall per 
unit area and including friction: 

dmf f d 3 p fi 



tot — r dr 



This has the same form as eq. ( |53|) for the driving force Fj r but with the replacement f? 9 — > 
Now, the ultra-relativistic case is particularly simple: to leading order in l/j w , the wall 
induces a sudden change in particle masses, mf + — > mf _, but leaves particle distribution 
functions as they were in the symmetric phase fi = fl\ (which are not the equilibrium ones 
in the broken phase). This allows the ^-integral in (|67|) to be performed and one obtains 



tot 
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where Am 2 = m 2 \~ and E iy+ in the momentum integral is also the one corresponding to the 
symmetric phase, as indicated. 

If F tot remains positive even for v — > 1, the system will enter the runaway regime. Besides 
the explicit form ( l68l) . this condition can be rephrased in terms of the free energies in the 
mean field approximation [27]. Notice that the mean field approximation implies 



(W)) MF = v^(o) + ^K 2 (0) 

i 

such that 

( V t)mf = ( V t)mf + J2 A 

i 

and the criterion for runaway behaviour reads 



dmf 



2 dV T 



dm 2 



(69) 



(70) 



< Ftot = {T + - F~) MF . (71) 

The results we just summarized can be interpreted as saying that the friction Ff r (which 

we expect to increase monotonically with v) saturates at a finite valu^l for v — > 1. A sketch 
for this behavior is given in Fig. [TTJ Noting that the driving force is given by 

F dr = 7- - JT+, (72) 

then the maximal value for the friction force reads ({F + ) MF = F + ) 

F T X = {F-)mf - r~ = ( V t')mf - V t (73) 



Using the convexity properties of the different contributions to Vt one can show that Fj^ ax 
is necessarily positive. 
If - T~) 

mf > 0' t nen Fdr < FJ^ ax , and the wall velocity will grow till the friction 
force equilibrates F& and a steady state with some terminal velocity is reached. In the 
opposite case with {J r+ — F~) MF < 0, one has F dr > F^ ax and the wall will keep accelerating 
without reaching a steady state, i.e. it will "run-away". The crucial assumptions here are 
that there are no hydrodynamic obstacles that prohibit the wall velocities to become highly 
relativistic in the first place, and that the mean free-path of the particles is larger than the 
wall thickness. 

It is instructive to rederive the results of ref. [27] we have just presented in the language 
of Kadanoff-Baym equations (see refs. [32] for an introduction to Kadanoff-Baym equations 
in Wigner space). In this formalism, the Wightman function Gf (for particle species i) is the 
relevant Green function that encodes the distribution function of the particles in the plasma. 
Under the assumption that the mean free-path of the particles that obtain a mass by the 
Higgs vacuum expectation value is much larger than the thickness of the bubble wall, one 
can neglect collisions in the wall. In this case the Kadanoff-Baym equation for the Wightman 
function Gf reads 

(p 2 - m\) e lo/2 G< = 0, (74) 



3 Up to a possible logarithmic increase with log 7™ (we thank Guy Moore for clarifications on this point). 
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< yv > 

Figure 11: A sketch of the friction force as a function of the wall velocity showing the saturation 
at (v) —> 1. The behavior for intermediate velocities is largely unknown. The arrows indicate two 
possible values for the driving force that would lead to steady or runaway bubble expansion, as 
indicated. 



where the diamond operator is defined as 

o = ^dx^p - d p d x . (75) 

In the semi-classical limit, the operator e i<>//2 can be expanded in gradients and the real/imaginary 
parts of the equations become, at first order, 

(p 2 -m*)G< = 0, (76) 
(p 2 -m*)oGf = 0. (77) 

The gradient expansion is justified in the present case because in the wall frame the particles 
have momenta of order jT, which is large compared to the inverse wall thickness. Using the 
ansatz 

Gf = AK5(p 2 -m*)f l (x» } p») , (78) 
this yields for the particle distribution function the equation 

1 



P • dx + ~(d c m?)\ ■ d p 



/,(^,^) = 0. (79) 



We introduced again a four-vector A M to parametrize the motion of the wall [that is, A M = 
(0, 0, 0, 1) in the frame moving with the wall] and ( = X-x. Notice that the first term in (1791) 
is the flow term of a Boltzmann equation while the second term represents the force from 
the wall acting on the plasma. In front of the wall the distribution function is given by the 
equilibrium one. E.g., for a bosonic degree of freedom 



Afro) 
exp[/3(w • p)] — 1 



fix = ^ a ,rL ; . (so) 
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where is the plasma velocity four-vector. Solving eq. ( I79p the distribution function behind 
the wall reads (with Am 2 as before) 

/, = -^£> (81) 
exp p \ I i — 1 



with 

Vti = u ■ p + u ■ A 
In the wall frame this yields 



A • p — sign(A ■ p) \ (A • p) 2 + Am 2 



(82) 



= IPo - iv J p 2 z + Am 2 . (83) 



This is in accord with the results of [27] and our previous discussion: in the wall frame the 
particles just cross the wall and change their momentum according to 

pi -> p 2 z - Ami (84) 

The contribution of particle species % to the pressure behind the wall is (in the wall frame 
and per degree of freedom) given by 

r rpplasma f j P 2 

0i ** ~ J (2nY Pz 1 

« 2 J ^~ 3 S(p 2 - m\ + ) p z ^p 2 z - Am 2 sign (p z ) 6(p 2 z - Am 2 ) f«, (85) 

where the last line is obtained by shifting the integration variable as in (j8%|) and noting 
that the measure dp z p z does not change. Hence, for the pressure difference produced by the 
plasma this gives 



SAT,, « Am? J 703<S(P 2 " m? + )/£ 



- ^Swh^> (86) 

which agrees with the formula obtained in ref. [27J. In terms of AT ZZ , the criterion for 
runaway solutions is very transparent and is simply 

e > AT ZZ = SiAT zz . (87) 

i 

In this case, only a part of the available vacuum energy is released into the plasma and the 
remaining energy is used to further accelerate the wall. 

This criterion for runaway behavior can be reformulated in terms of parameters of the 
phase transition in case it is rather strong. For particles that are light in both phases, the 
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contribution to the pressure is similar in both phases. Particles that are heavy in both 
phases do not contribute much to either pressure difference or free energy. Hence, mostly 
the particles that become heavy during the phase transition produce a pressure difference 
along the wall. Using this, one obtains for an accelerated wall the criterion 

light— theavy light— theavy 

with Cj = 1 (1/2) for bosons (fermions), |A^| are the corresponding numbers of degrees 
of freedom, are the coupling strengths to the Higgs boson and (0) the Higgs vacuum 
expectation value in the broken phase. Using the relation e = ^(oatT^) we get that a 
runaway wall is in principle possible for 

30 ^2light^heavy Ci 1^1 Vi /onX 

a N >a 00 = — [ — \ y —- , 89 

\T N J 22 ligM Ci\Ni\ 

with d i = 1 (7/8) for bosons (fermions). This equation serves as the definition of a^. In 
extensions of the SM, eventually more particles contribute to the pressure difference, but 
typically not many new light particles are in thermal equilibrium. One hence can deduce 
that for 

a N > 1.5 x 10- 2 fMj , (90) 

runaway walls are possible depending on the details of the model. It is interesting to note 
that models which lead to sizable gravitational wave production typically satisfy the runaway 
condition and this should be taken into account when calculating the GW signal. 

Next, we make contact with the case of a constant wall velocity, discussed in the last 
section. 



7 Energy budget of first-order phase transitions 

The analysis in the last section assumed that the system was time-independent in the wall 
frame, which leads to the fact that the Higgs field only contributes a pressure component 
from the vacuum energy to the energy momentum tensor 

T*\t = -e, I&£ = 0. (91) 

In a static system these contributions have to be compensated by the plasma and this requires 

rpplasma |+ ^ rpplasma |+ g (92) 

Such relations lead to the matching conditions used as boundary conditions in the hydro- 
dynamic analysis of the plasma. In the case of a highly relativistic plasma, these boundary 
conditions can be derived explicitly from the particle distribution functions. Following the 
same calculation as in the previous section, for j^ asma one finds 

rpplasma |+ rpplasma |+ g (93) 



28 



T| = 0.2 T| = 1.0 

deflagration detonation run-away deflagration detonation run-away 




0.1 1 10 0.1 1 10 

Figure 12: The energy budget for r\ = 0.2 and rj = 1.0. The different contributions (from top 
to bottom) are thermal energy, bulk fluid motion and energy in the Higgs field. The last two 
components can potentially produce anisotropic stress in the plasma and subsequently gravity 
waves. 



Hence, in the runaway case, with > the solutions for the fluid motion are identical 
to the ones with = «oo, according to the distribution functions determined close to the 
wall. At the same time the Higgs field cannot be time-independent anymore and energy 
momentum conservation implies that the remaining energy is used to accelerate the wall. 

We observed in section 4 that, in the limit of large wall velocities, the efficiency factor 
does not depend on the wall velocity but is given by (j4"2"j) . This means that, in the runaway 
case, 

(runaway) . (94) 



0.73 + 0.083^0^ + a 

In summary, in the runaway regime and for given a^, a portion of the initial produces 
bulk motion with efficiency k^, as given by eq. ( I9~4"j) . while the remaining portion, ot^ — a^, 
is transformed directly into kinetic/gradient energy of the Higgs field with efficiency k = 
1. These two components can potentially produce anisotropic stress in the plasma and 
subsequently gravity waves while the thermal energy in the plasma can not. Figure [12 
shows the energy budget of the phase transition for two choices of the friction coefficient rj 
as a function of in different regimes of bubble expansion. 



8 Summary 

The bubble wall velocity £ w in first-order phase transitions is a key quantity entering the 
calculation of the baryon asymmetry in electroweak baryogenesis and its derivation has been 
discussed extensively in the literature. However, it has been treated in detail only in specific 
models (corresponding to weak first-order phase transitions) and a general account of the 
problem was lacking. In this work, we attempted to gather all the important information 
in a self-consistent manner and in a model-independent approach. We presented a unified 
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description of all the different regimes characterizing bubble growth and stressed how they 
are connected. 

In the last few years, there has been a significant effort towards working out the relic 
gravity wave background generated by bubble collisions. Production of gravity waves during 
a first-order phase transition is due to bulk motions in the plasma and magnetic fields 
generated by these fluid motions. The resulting gravity wave spectrum scales as a large 
power of the typical fluid velocity, roughly as v A . The fluid velocity profile in the vicinity 
of the bubble wall is relevant and it is therefore important to estimate the fraction k of 
the vacuum energy density e liberated during the phase transition which goes into the fluid 
motions. 

So, £ w and k are the two salient quantities that one would like to predict for any given 
particle physics model leading to a first-order phase transition. The most straightforward 
quantity to compute in a given model is a at, the ratio of the vacuum energy density e to the 
radiation energy density at the nucleation temperature. One would be interested in knowing 
in which regime (detonation, deflagration, hybrids, runaway) the given model is expected to 
fall in, without having to deal with the intricate Boltzmann equations. There is a relation 
which has been used extensively in the literature for £ w as a function of a^. Such relation is 
valid for supersonic walls and is based on the assumption of Jouguet detonations. Assuming 
Jouguet detonations amounts to setting the value of the fluid velocity at the inner boundary 
of the velocity profile in the wall frame to the speed of sound. There is no justification for 
this choice although it has been used in the literature for simplicity. Setting arbitrarily some 
boundary condition for the velocity leads to inconsistencies and corresponds to ignoring the 
constraints from the equation of motion of the Higgs field. By dealing explicitly with that 
Higgs equation, our study goes beyond this assumption. 

The difficulty is that the problem under consideration is not fully determined by a^. As 
we elucidated, for a given a at, there are many possible wall velocities and a large range of k 
values^, as seen in Fig. |EJ To fix the solution, one needs to compute the friction term that 
restrains the bubble expansion. There are two approaches to this problem. The rigorous 
(but tedious) treatment requires solving coupled Boltzmann equations. A simpler approach 
is to take the friction as a parameter independent from the bubble wall velocity and solve 
the equation of motion for the Higgs field by using some phenomenological description for 
the friction. This is the approach that we took in this paper. It has been adopted in some 
numerical studies [T ^[T6T - [T9] . although the modeling of the friction was not appropriate in 
the limit of large velocities. We now use a more realistic modeling which allows us to describe 
the runaway regime in addition to the steady state regime, eq. ( H9|) . 

One of our final results has been to present, in Fig. [TOl model-independent contours for 
£ w and k in the r\ — plane, where r\ and are the two crucial physical (dimensionless) 
parameters respectively characterizing the strength of the phase transition and the amount 
of friction. In concrete models, oln and r\ are not independent parameters. However, we 
believe that the information contained in these plots is useful as it enables to investigate the 
parametric dependence of the wall velocity and the k factor in different regimes and gather 
all the important physics in a single comprehensive figure. The contour plots are not fully 
model-independent as we have fixed the relative change in the number of relativistic degrees 

7 We have also corrected a missing factor in the definition of K which is found in the literature. 
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of freedom between the two phases but they serve the purpose of describing the qualitative 
behavior of £ w and k. From our analysis, we obtain that only a small band in the (77, ajv) 
plane leads to pure detonations. We also find that in a large region of this plane, hybrids 
coexist with runaway solutions, meaning that our time-independent analysis is not sufficient 
to fix the solution. According to numerical studies, hybrids are mostly unstable [TH]. We 
would therefore conclude that the most likely solutions are either deflagrations or runaway 
solutions. We have clarified the condition determining the runaway regime, see eq. (|89|) . 
and found that for values ~ 0(1) common in the literature, a runaway regime is quite 
likely. Note however than numerical studies have not been carried out in the ultra relativistic 
regime. 

For deflagrations and hybrids, a significant fraction of the available energy goes into bulk 
fluid motions. On the other hand, for a non-steady state solution, a large fraction of the 
energy goes into accelerating the wall and little goes into bulk fluid motions, as seen in 
Fig. [T2j Therefore, and perhaps counterintuitively, we find that, in a very strong first order 
phase transition, the contribution to the gravity wave spectrum from turbulent fluid flows 
is probably subdominant. While the overall size of the gravity wave signal is essentially 
controlled by the amount of vacuum energy released, i.e. by a^, its detailed spectrum 
will depend on how the energy is distributed among the different components (wall versus 
plasma). We have therefore provided the relation between the fluid velocity and the bubble 
wall velocity in Fig. [7J In addition, we estimated the thickness of the plasma shell near 
the bubble wall where the kinetic energy in the plasma is concentrated. In all numerical 
calculations of the GW background [HEHE], it has been assumed for simplicity that all the 
energy was concentrated on the bubble surface and finite thickness effects were ignored (they 
were instead considered in the analytical approach of [9]). This is a reasonable approximation 
for detonations, although we show that even for supersonic wall velocities, the thickness of 
the plasma shell can reach ~ 20% of the bubble size, see Fig. E0 

Finally, the nature of the electroweak phase transition is unknown and it will take some 
time before we can determine whether electroweak symmetry breaking is purely Standard 
Model-like or there are large deviations in the Higgs sector which could have led to a first- 
order phase transition. Although our analysis was implicitly motivated by the electroweak 
phase transition, it could be applicable to other first-order phase transitions in the early 
universe. 
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A Numerical fits to the efficiency coefficients 

In this section we provide fits to the numerical results of section HI These fits facilitate the 
functions k(£ w , ajy) and a + (£ w , a^) without solving the flow equations and with a precision 
better that 15% in the region 10~ 3 < on < 10. 

In order to fit the function k(£ w , ckjv), we split the parameter space into three regions 
and provide approximations for the four boundary cases and three families of functions that 
interpolate in-between: For small wall velocities one obtains (£ w <C c s ) 

K A — ^w b ~r^a n nQ7 r^^i ' (^5) 

1.6b — 0.037 waN + o>n 



For the transition from subsonic to supersonic deflagrations (£ w = c 

2/5 

0.017+ (0.997 + tt7v) 2/5 ' 
For Jouguet detonations (£ w = as stated in eq. (BE 



«b - 777^-7^—^ ■ 06) 



Kn — = and tj = . 97 

0.135 + ^0.98 + ^ l + a N v ' 



And finally for very large wall velocity, (£ w — >■ 1) as stated in eq. (112]) 

«u — • 98) 

0.73 + 0.083v7i^ + ajv 

For subsonic deflagrations a good fit to the numerical results is provided by 

11/5 

[C s ~i w ! )K B +£ w Cs' K A 

and for detonations by 

(6 - i) 3 ey v /2 «c« g 

[(6 " l) 3 - - l) 3 ]d /2 ^c + - 1) 3 «:d 



> o) - c " ... • (ioo) 



The numerical result for the hybrid (supersonic deflagration) region is well described by a 
cubic polynomial. As boundary conditions, one best uses the two values of k and the first 
derivative of k at £ w = c s . Notice that the derivative of k in £ w is not continuous at the 
point £j. The derivative at £ w = c s is approximately given by 

5k -0.9 log ^ (101) 
1 + \/a/v 



This differs from the derivative one would obtain from the fit in the region £ w < c s , but 
mostly for values a > 1, where no solutions exist for £ w < c s . The expression for supersonic 
deflagrations then reads 

«(c a < £ w < O) - k b + - c s )Sk + — -^[k c - «b - (6 - c a )<*«] . (102) 

(O - c s ) d 
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The fits of k compared to the numerical results are given in Fig. [13j The relative errors 
never exceed 15%. 

Another useful function is the maximal a/v that can be realized for fixed £ w . This function 
is approximately given by 

a nu»„l( 1 _ ^-13/10. (103) 

Finally, the connection between and a + might be needed in phenomenological studies. 
Like in the case for k, we provide approximations or analytic expressions in limiting cases 
and interpolate in between. For velocities larger than the Jouguet velocity £j, a + and 
coincide. For < 1/3 all velocities can be realized and in this case the limit v w — > also 
yields a + — > a^. Otherwise there is a minimal £ w which is, according to eq. f)103j) . given by 

a+ = 1/3 for £ w = 1 - (3a N y w/13 if a N > -. (104) 

For relatively small a wall velocity of the speed of sound can be realized leading to a fit 
a + {c s ) ~ ajv(0.329 - 0.0793 log a N + 0.0116 log 2 a N + 0.00159 log 3 a N ) 

if a N < -(1 - c s )- 13/10 ~ 1.02. (105) 
3 

Finally, the transition from supersonic deflagrations to detonations can be determined ana- 
lytically leading to 

3 — 3£ 2 

a+(6) = 7 for a11 aN ' ( 106 ) 

For (y 7v > 1/3, the function a + (aN,£, w ) is well approximated by a polynomial in £ w that 
contains the two/three data points just given for the case ^ 1-02. For ajy < 1/3, the 
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function a + (ajy, £ w ) is approximately linear in the region £ w e [c s , £j] and goes slowly to ajv 
for £ w — >■ 0. For velocities below c s , a reasonable approximation is of the form 



and matching of this function and its derivative to the linear regime yields 

( c s - - a+(c s )]a N - 2£*(1 - Cj/c s )[a + (c s ) - a7v]a+(c s 



(107) 



a: i 



( C 2 - - a+ic)] - 2C(1 - Ucs)[a+(c s 



a N 



(108) 



The comparison between fit and numerical results is given in Fig. [TH The relative error 
never exceeds 5% and is worst for a at ~ 0.3. 




Figure 14: The ratio a+/aN as a function of the wall velocity £ w for fixed ck/v and the fit (dashed 
lines) described in the text. 
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